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We present analytical and numerical evidence for the validity of an effective S^s 



approach 



to the description of random field generation in S > 1, and especially in an 5 = 1, dipolar spin 
glass models with strong uniaxial Ising anisotropy and subject to weak external magnetic field Bx 
transverse to the Ising direction. Explicitely Bx— dependent random fields are shown to naturally 
emerge in the effective low-energy description of a microscopic S = 1 toy model. We discuss our 
results in relation to recent theoretical studies pertaining to the topic of B^;— induced random fields 
in the LiHoa;Yi_a;F4 magnetic materials with the Ho^^ Ising moments subject to a transverse field. 
We show that the S^s — ^ approach is able to capture both the qualitative and quantitative aspects 
of the physics at small B^ , giving results that agree with those obtained using conventional second 
order perturbation theory. 



I. INTRODUCTION 



In condensed matter physics systems with strongly in- 
teracting quantum mechanical degrees of freedom, it is 
often a challenge to explain physical phenomena from a 
truly first principle atomistic point of view. In systems 
where there are high energy scales well separated from 
a low-energy sector, effective low-energy theories offer 
the advantage of a reformulation of the problem with 
an exponentially smaller Hilbert space. A well known 
and topical example where such an approach is used is 
in the derivation of an effective spin-only model starting 
from a Hubbard model describing electrons hopping on 
a lattice. It is commonly accepted that the low-energy 
magnetic excitations of a Hubbard model with a large 
Coulomb repulsion U are easier to investigate within an 
effective spin HamiltonianFii^i^ Generally speaking, the 
only requirement to be able to derive an effective model 
is to have a small parameter, which is t/U in the pre- 
vious example, where t is the nearest-neighbor hopping 
constant. 

In many magnetic materials, the ground state degen- 
eracy of the otherwise free magnetic ions can be partially 
lifted by electrostatic and covalent interactions due to the 
surrounding atoms — the so called crystal field effect. In 
a number of situations, the energy scales associated with 
the spin-spin interactions are much smaller than the en- 
ergy gap between the single-ion ground state and the 
excited crystal field states. In such cases, one can, as a 
first approximation, often neglect the high energy states 
and reduce the relevant Hilbert space to a much smaller 
subspace of low energy states. In this paper, we dis- 
cuss the quantitative validity of an effective low-energy 
theory description of a model inspired by the phenom- 
ena displayed by the disordered LiHo3;Yi_2;F4 magnetic 
material when subject to an external magnetic field 
applied perpendicular to the Ising direction of the Ho'^''' 
magnetic moments. 



The LiHo2;Yi_2:F4 magnetic material exhibits many 
interesting magnetic behaviors i^i^ii'^'^iiiS The magnetic 
properties of LiHo2;Yi_a;F4 are due to the Ho^+ ions. 
The single- ion ground state of Ho^^^ is a doublet, 
while the first excited state is at ^ 11 K above the 
ground stateJ^ii^ The most relevant interactions be- 
tween the magnetic Ho'^+ ions are magnetic dipole-dipole 
interactions]^ Since the maximum strength of the dipo- 
lar interactions is for nearest neighbor separation and is 
approximatly 0.31 K, collective behavior in this mate- 
rial occurs at temperatures less than 0(1 K) where only 
the ground doublet is significantly thermally populated. 
Consequently, the cooperative phenomena and the low 
temperature properties of this material in zero applied 
magnetic field should be well captured by an effective 
model with spin-i degrees of freedomJ^iii For example, 
in zero applied magnetic field, the system can be recast as 
a diluted dipolar Ising model with the low-temperature 
phase being either a ferromagnet or a spin glass depend- 
ing on the concentration x of magnetic ions4iii^ On the 
other hand, for a; = 1 and with a magnetic field ap- 
plied perpendicular to the crystallographic Ising c— axis 
direction, LiHoF4 has been advocated as one of the rare 
physical realization^ of the transverse field Ising model 
(TFIM),i^iiiii^ii^i^ Yet, it is only relatively recently that 
a somewhat rigorous justification of a TFIM description 
of LiHoF4 in nonzero B^ has been put forwardJ^ How- 
ever, over the past twenty years, and until very recently, 
several experimental studies had found the behavior of 
LiHo3;Yi_3;F4 (x < 1, Bx =/= 0) paradoxical, as we now 
discuss. 

One may have naively expected that the application 
of a transverse magnetic field in LiHoa;Yi_2;F4 would 
allow to explore the physics of the TFIM in either a 
diluted ferromagnet or a spin glass, depending on the 
concentration x. However, the situation for x < 1 is 
quite a bit more complicated.— For example, for B^ = 
0, LiHoo.i67Yo.833F4 displays a conventional spin glass 
phase transition^i2^ with a nonlinear magnetic suscepti- 
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bility, X3 , diverging at the spin glass transition tempera- 
ture, Tg, as XaiT) cx (T — Tg)~^ as in ordinary spin glass 
materials.— However, as is increased from zero, Xz{T) 
becomes steadily less singular, and there appears to be no 
Ba;— induced quantum critical phase transition between a 
paramagnet and a spin glass statej^ This puzzling exper- 
imental behavior had been tentatively interpreted as due 
to a 1^' order transition near the T = quantum phase 
transitioni^i^ However, very recent and independent the- 
oretical investigations2^i2£i^ have instead proposed that 
the microscopic origin of the "quenching" of the param- 
agnetic to spin glass transition as B^ is turned on is due 
to the generation of random fields that destroy the spin 
glass phase. 

The authors of Ref. [1^ used an effective S'eff = f/2 
theory, very similar to the one developped for pure 
LiHoF^i^ to expose how random fields develop in a mi- 
croscopic model of LiHoa;Yi_3;F4 in nonzero Bx- In par- 
ticular, Ref. [2^ showed how the nonlinear susceptibil- 
ity X3 becomes progressively less singular as Bx is in- 
creased. Also motivated by the phenomena displayed by 
LiHo2;Yi_a;F4, Schcchtcr and coUaborator a^^i^^'^^ also re- 
cently investigated in a series of papers the general phe- 
nomenology of induced random fields in LiHo2^Yi_a;F4. 
To do so, they considered in Refs. [25lf27j an easy-axis 
spin-S' {S > 1) dipolar spin glass toy model Hamiltonian, 
H, in presence of a nonzero Bx- By using second order 
perturbation theory, invoking the scaling droplet picture 
of Fisher and Huse for spin glasses and using an Imry- 
Ma type argument;^ Schechter et al^'^ calculated the 
finite energy SE required to flip the spins within a spin 
glass droplet, finding a limit on how large the spin glass 
correlation length ^ can grow to as the system is cooled 
from the paramagnetic phase. The behavior of the sys- 
tem, and the corresponding SE, is found to be analogous 
to that of a spin glass in a random magnetic field which, 
according to the droplet model, does not show a spin g lass 
transition in nonzero field.— As a result, Refs. [25ll27| ar- 
gue that no spin glass transition can occur in a dipolar 
spin glass where random off-diagonal dipolar interactions 
and an applied transverse magnetic field are simultane- 
ously at play. 

On one hand, the results of both Refs. [25ll27| and 
Ref. [2^ derive from the notion that, the applied trans- 
verse field generates, through the off-diagonal part of the 
dipolar interactions, which couple the Ising z component 
with the perpendicular x and y components, some effec- 
tive random fields. However, it has so far not been clar- 
ified to what extent the random fields are quantitatively 
equivalent or only qualititatively related in those two sets 
of works. In their studies, the authors of Refs. (25ll27t ar- 
gued, correctly, that considerations of a model with large 
spin (5 > 1) is crucial to understand the weak field re- 
sponse of the spin glass phase in either their toy model H 
or in LiHoa;Yi_a;F4. Exact diagonalization results of an 
5 = 1 dipolar spin glass model with easy-axis anisotropy 
provided further quantitative support to the theoretical 
arguments as to the scaling behavior of 5E with both Bx 



and the number of spins in the systemi^^i^ At the same 
time, their results from similar calculations^! for an ef- 
fective anisotropic spin-i dipolar Ising model in a trans- 
verse field, but with the off-diagonal dipolar interactions 
rescaled compared to the longitudinal Ising coupling;^ 
did not conform with those obtained for the "bare" (high- 
energy) anisotropic S = 1 model.— i^I Partially on the 
basis of those results, and seemingly confirming a pre- 
vious argumeni*^^, Ref. [l^j, concludes that an effective 
spin-i model, such as that used in Ref. [2^, is not suffi- 
cient to capture the physics in the small Bx regime com- 
pared to the "bare" microscopic (large-spin) anisotropic 
dipolar spin glass model H. The question of the use- 
fulness of an effective spin-^ model to describe random 
field phenomena in the dilute ferromagnetic regime of 
LiHoajYi-ajF^ i^^i"^^ has also been recently raised<^ 

Considering a perspective beyond the specific problem- 
atic of LiHoj;Yi_2;F4, one could interpret the conclusion 
of Refs. ^25ii27li28t regarding the inadequacies of an effec- 
tive spin-i model to describe LiHo3:Yi_^F4 in Bx 7^ 
as a counter example of the precise quantitative useful- 
ness of effective low-energy theories for quantum iV-body 
systems. It is therefore useful to investigate with some 
scrutiny the mathematical justification for an effective 
spin-i model for LiHoa:Yi_a:F4 with Bx ^ 0. This is the 
purpose of the present paper. More specifically, the ques- 
tion that we ask here is: to what extent arc the cxplicitely 
manifest random fields derived in an effective low-energy 
theory, such as in Ref. [2^, related to the random field 
like effects at play in perturbation theories, such as used 
in Refs. [25l|27| ? Below we show, via a derivation of an 
effective low-energy S'eff = ^ Hamiltonian for anisotropic 
dipolar glasses, that effective random longitudinal fields 
emerge naturally in the SeS = \ model. On the basis 
of analytical calculations and exact diagonalizations, we 
highlight the fact that an SeS — \ Hamiltonian properly 
derived from an 5 = 1 high-e nerg y toy model 7i, such as 
the one proposed in Refs. (25lf27f (see Eq. (1) in Section 
ITT|) . is a quantitatively valid and controlled approach to 
this problem. 

The paper is organized as follows. We first discuss 
in Section |TT] an anisotropic spin-5 dipolar Hamiltonian 
as a simplified model displaying the key physics of the 
LiHOj.Y1_j.F4 material in a transverse field and show in 
Section IIIII how to derive from it an effective S'eff = 1^ 
Hamiltonian to lowest order. We present in Section IIVI 
results from exact diagonalization calculations that com- 
pare the S = 1 and the Seff = \ models and which di- 
rectly confirm the quantitative validity of the effective 
Hamiltonian approach. Section fVl concludes the paper. 



II. ANISOTROPIC SPIN HAMILTONIAN 

The Ho'^+ ion is characterized by a very large hyper- 
fine interaction between the electronic and nuclear mo- 
ments and the effects of this strong interaction plays an 
important role in a number of Ho^+— based magnetic 
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materials4ii^iiii2^i^ In particular, in LiHoF4, it leads 
to a significant increase of the zero temperature criti- 
cal transverse field for the dipolar ferromagnet to quan- 
tum paramagnet transition. ^'^i^^ It also plays an impor- 
tant role in setting the relevant critical transverse mag- 
netic field scale in the dilute LiHoxYi_i:F4.— In this pa- 
per, however, we arc specifically interested in the gen- 
eral phenomenology of random fields along the Ising spin 
directions generated by small applied transverse field 
rather than obtaining a precise quantitative description 
of LiHo2:Yi_a,F4. In this specific context, we therefore 
neglect the role of hyperfine interactions. Also neglecting 
the hyperfine interactions, Schechter et al^^^ proposed 
a generic anisotropic spin-S* toy model Hamiltonian with 
long-range dipolar interactions 
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This Hamiltonian is a simplified model that preserves the 
basic characteristics of the proposed microscopic Hamil- 
tonianHi^ for LiHoa;Yi_j;F4. In the absence of an ex- 
ternal field, individual IIo^+ spins have an Ising like 
ground state doublet with a large energy gap between 
the next excited state and the ground doublet. Also, 
for S* = 1, the excited state of model in Eq. ([1]) is a 
singlet, as for 110"^+ in LiHo3;Yi-aF4 i^^'^^ Here, i,j are 
the positions of the randomly positioned magnetic mo- 
ments. Vf^'^ denotes the random long-range dipolar in- 
teraction between the spins, where T^J^ stands for the 
Ising interaction and V/j^ stands for the off-diagonal in- 
teraction {V^" = V^'^ for dipolar interactions). D > 
is the anisotropy constant mimicking the crystal field. 
For Bx = 0, the ground state (GS) of a single spin is 
doubly degenerate with ~ ±5 . The correspond- 
ing states of the doublet are denoted \S) and | — S"). 
The first excited states have = zt{S — 1) and energy 
Hq = (25'— 1)13, with the corresponding states denoted as 
I ± (5 — 1)). Ignoring momentarily the V^'^'^ interactions, 
the Zeeman term, —Bx Sf, lifts the GS degeneracy of 
the I ± S) ground doublet, resulting in two new lowest 
energy states, \a{Bx)) and \P{Bx)), with corresponding 
energies Ea{Bx) and Efj{Bx), and with an energy gap 



^{Bx) — Ea{Bx) — Ep{Bx 
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between them. For B^ <^ flo, to leading order in pertur- 
bation theory, the gap A{Bx) is proportional to {Bxf'^ ^ 
Invoking the spin glass droplet scaling picture of Fisher 
and Huse^^ and using an Imry-Mai^ type argument, one 
can calculate the energy required to flip a spin glass 
droplet of size L containing N ^ L"^ spins, with d the 
number of space dimensions (here d = 3). This en- 
ergy cost is due to the perturbative quantum T-l± = 
— Y^i^j Vfj^S^Sj — Bx which term does not com- 

mute with the the classical = -DY.^[[Slf - S'^] - 



^Yi^jVfj^SfSj term. Considering first only T-i\\, and 
taking the droplet picture of only two distinct ground 
states,— |$5) and |$s) denote the collective (doubly- 
degenerate) Ising spin glass ground states of the sys- 
tem. These two ground states are related by the global 
Sf — > ^Sf symmetry, where each spin is either in its 
\+S) state or its \~S) state. As discussed in Refs. 25ll27j . 
nonzero Ti.± lifts the ground state degeneracy, as we now 
review in order to make contact with the results pre- 
sented below in Sections IIIII and IIVI 

The lowest energy excited states (above the otherwise 
two degenerate |<i>5) and \^s) ground states) are 
and |0(^5_i)) states, in which the fc'th spin has its 
quantum value changed from +S to -|-(S' — 1) or from 
—5* to —(5—1). Using standard second order degenerate 
perturbation theory^ and considering only excitations 
to the (intermediate excited) \(f>(^s-i)) l<^(s_i)) states, 
the fluctuation-induced energy difference between |<I>s) 
and \^s) is 



SE = 



-H. 



4|^*s.$.J' (3) 



where 
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^0 V 
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<^s\Hi_\ 



*fs-i) 



where we have taken the ground state energy to be zero. 
Since ($s \Hx_\ 4>\s-i)) = <^fs-i)) = 0' l^^ve 



■^*s,$s = Subtracting H^^^^^ from ff*s,<[.s, only the 
odd terms in Bx remain, with the even terms in Bx can- 
celling each other out. Finally, to lowest order in Bx, we 
get 



(4) 



Taking the largest V^^^' with a typical value Vj^ , the typ- 
ical energy gained by flipping a droplet oi N ^ spins 
is, to leading order in Bx, 



S'BxV^^N 

no 



(5) 



indicating that the total energy gain increases with Bx 
linearly to leading order, as first found in Refs. [25ll27| . 

This decrease in energy is to be compared with the en- 
ergy cost due to the formation of a spin glass droplet 
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This energy cost scales with the hnear size L of the 
droplet, L = N^/^, as « S^Vi^L^", where V|| is the typi- 
cal value of the largest Vfj^, which one typically expects 
to be of the same order as V±. Comparing the energy 



gain (\SE\) of Eq^ 
formation, Refs. [25l |' 



with the energy cost for droplet 
I find a finite correlation length ^, 



identified with L, which, for small Bx, scales as 



(6) 



Based on an argument by Fisher and Huse^ 9^ < 
{d — l)/2, or 9d < 3/2 here. Hence, turning on leads 
to a reduction of the correlation length £^{Bx), inhibiting 
its divergence as occurs when B^ = 0. In other words, 
the presence of the applied transverse B^ leads, via the 
presence of the off-diagonal spin-spin interactions, to 
a destruction of the spin glass phase with a typical spin 
glass correlation length ^ decreasing as B^ increases. As 
argued in Refs. 25ll27| . this is the mechanism via which 
the non-linear magnetic susceptibility X3 no longer di- 
verges in LiHoxYi_a;F4 as B^ is increased from zeroi^i^^ 



III. EFFECTIVE SPIN-i DESCRIPTION 



In the previous section we reviewed the arguments of 
Refs. [25II27I I which lead to the key result of Eq. (g]). We 
now proceed to show that a reformulation of the micro- 
scopic spin Hamiltonian, Eq. ([1]) , in terms of an effective 
Scff ~ \ model, leads identically to Eq. ^ in the limit 
of small B.J./D. 

Firstly, we focus on a situation where the temperature 
considered is low compared to Hq, and project the spin S 
operators onto the two-dimensional subspace formed by 
the two lowest energy eigenstates, \a[Bx)) and \ j3[Bx))- 
Following Refs. [Ullal, we define an Ising subspace, 1 1) 
and I J,), by performing a rotation 

IT) = -^(|a>+cxp(z0)|/3)) 

U> = ^(|a)-cxp(z0)|/3)). (7) 

The phase is chosen such that the matrix elements of 
the operator within the new (Ising) subspace are real 
and diagonal. In this case, we can define Sf = Czz^^l- 
This allows us to recast H. in Eq. ([T]) in terms of an 
effective spin— ^ Hamiltonian, Tioff, that involves the 
Pauli matrices.— In this projected subspace, a transverse 
field r = \l\(Bx) acts on the effective uf spin. The 
projected {fi = x, y, z) operator may be written as: 



(8) 



The C^j/ and A dependence on B^ can be obtained by 
exact diagonalizationi^i^ of the non-interacting part of 
n (i.e. Vl^" = 0) in Eq. ©. 
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FIG. 1: (Color Online) Evolution of A, Czz, CxO, and Cxx as 
a function of the external transverse field Bx for S = 1. 



For zero transverse field, B^ ~ 0, the only nonzero Cfn, 
coefficient is Czz{Q) = S, giving a ^'classical" (effective) 
low-energy dipolar Ising model 



Hising = ^^^^ X! ^Ki'^'^i'^j 



(9) 



Turning on B^, the coefficients Cxo and Cxx increase 
with Bx, while C^z shows a slight decrease with increas- 
ing Bx, as shown in Fig. [TJ Thus, by substituting S'f 
with CzziBx)<y! and Sf with C(B,)af + C,o(S,)l in 
Eq. Il]), the effective spin-i Hamiltonian is 



Tic 



1 



Clz{Bx 



-Czz{Bx)\Cxx{Bx)Y.Vt^ata] 



(10) 



- Cxo[Bx)Y,V,fat I - \a{Bx)Y,cj^ . 

ij^j ) i 

As can be seen, the projection of the Vfj^S^Sj term in 
Eq. ([T]) results in an induced random bilinear coupling, 
cx o'^a^, and a longitudinal random field interaction, oc 
tjf , for Bx 7^ 0. For low enough transverse field Bx, the 
Ising dipolar interaction (c>c Vfj^) is the dominant term. 

Having derived the effective Hamiltonian, we now re- 
peat the calculation of 6E within this effective ^off = ^ 
framework by again bringing in the spin glass droplet 
picture^S. For Bx = 0, we denote the ground state 
of the SeS = \ system where is a specific realiza- 
tion of the t and [ (effective) Ising spins configurationi^ 
For Bx = 0, because of time reversal symmetry, the time 
reversed state \ip), which is obtained by ffipping all the 
spins of 1^/;), is a ground state of the system as well, giv- 
ing a ground state doublet in the "effective spin" droplet 
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picture. Carrying on a similar discussion as in the pre- 
vious section and as in Refs. [1^113, at low enough 
within a droplet picture, the symmetry is broken due the 
presence of the induced random fields in Eq. pU]) . The 
energy cost to flip the spins over a droplet is, 

5E = (^|i/off|V^> - (V'|i?cfr|V'> 
which, to lowest order in B^, gives 



6E « 2C,,CMJ2^^T 



(11) 



Although we have an exact analytical expression for the 
Cfj,i, coefficients as a function of B^ (which is available 
for S < 3/2), in order to compare with Eq. (|4]) above 
and with Refs. p5ll27l |. we consider the Bx dependence 
of the C^i/ to leading order in Bx/D. Using standard 
degenerate perturbation theory, for 5 > 1^^ the 1 1) and 
1 1) defined in Eq. ([7]) are, up to second order in Bx, given 
by 



IT) = 1- 



40§ 



s 



\S) 

\-s) + 



Bx /<5 
f^o V 2 

Bx 



S-l) 

I- 



recalling that Hq 
via which the C^i 



= (2s - 1)D. 
are obtained, 



(12) 

Returning to Eq. ([8|), 
e.g., C., = i((T|51T 
we use 



) - a \S^\i)) and Cxo = i((T 1^1 T) + (i 1^1 D), 
Eq. (dSl) to findC,, « 5(1 -||), Cxo « SBx/no, Cxx oc 

(S,)2S-i (Cxx « Bx/no for S = 1), while A cx (Bx)^^ M 
Substituting those Bx dependencies back in Eq. (fTTj) , the 
dependence of the energy cost SE is, to lowest order in 
Bx, 

'^^E^^r(V'Ki^) ■ (13) 



5E 



2S^ 



As we can see, the energy cost obtained in the Scs 
picture is identical to the energy cost given by Eq. ([4] 
obtained via second order perturbation theory and pre- 
viously reported in Refs. [25;27]. Thus, Eq. ([T3|) leads to 
the same RMS energy cost for flipping a droplet, given by 
Eq. (O, and the same Bx dependence of the spin glass 
correlation length ^ in Eq. ([6]). Hence, we have shown 
that a formally derived effective S'cff = 1/2 Hamiltonian 
does capture quantitatively the low energy physics of the 
full S Hamiltonian at low transverse fields. While the 
argument above was constructed for the toy model of 
Eq. ([T]), one could proceed identically for the full blown 
microscopic Hamiltonian of LiHo2-Yi_a;F4. Indeed, this 
is what is the underlying program of Ref . [26l | . 



IV. NUMERICAL RESULTS 



In the same spirit as Ref. 25ll27| . in order to inves- 
tigate to what extent our proposed low energy effective 



spin-i model is a good description of the full anisotropic 
Hamiltonian ([1]), and to determine the range of trans- 
verse field over which the above analytical small Bx field 
results is valid, we have performed numerical calculations 
to backup our perturbative approach. In this section we 
present results from exact diagonalizations on finite-size 
clusters with open boundary conditions^. In order to 
compare the present approach with the previous inves- 
tigations done by Schechter et al.r^^ we work at the 
same constant dipole concentration x = 18.75%. 

LiHoF4 is a compound with space-group C^^ {Mi /a) 
with lattice parameters a = b = 5.175A, c = 10.754, 
and has 4 holmium ions per unit cell positioned at 
(0,0,1/2), (0,1/2,3/4), (1/2,1/2,0) and (1/2,0,1/4)4^ 
For LiHoj;Yi_a;F4, a dilution of a; = 18.75% is realized by 
distributing randomly N magnetic moments (holmium, 
Ho'^^ ions) in a sample of possible sites. We have 

chosen samples of size (2a, 26, c x y ), where TV is a multi- 
ple of 3. Thus, changing the number N of magnetic ions 
means changing the size of the sample in the z-direction 
in order to keep a constant dilution. 

In Eq. ([T]), the dipolar interaction is written as V^^ , 
which takes, with the negative coefficient convention used 
in Eq. ([T]), the explicit form: 



Q/3 



3q0 



(14) 



where r.y is the distance between the ions at positions i 
and j, and a,[3 ~ x,y,z. The dipolar interaction V^^ 

is of the order ^ w 4.49 x lO"'^ K, whereas the on-site 
anisotropy is taken as D = 10 K. In the following, we in- 
vestigate the behavior of the gap SE between the ground- 
state and first excited-state as a function of the applied 
transverse field Bx. Since we are mainly interested in 
checking the relations ([U, ^ and (fT3|) . we present our 



results in terms of renormalized parameters 



SE 



)• 



To perform a first check of the validity of our approach, 
we choose a small cluster with a fixed random distribu- 
tion of A'^ = 9 spins and compute the renormalized gap 
5E/{dVN) for both models {i.e. S ^ 1, Eq. ^ and 
S'cff = |, Eq. pO)) ) as a function of the reduced transverse 
magnetic field Bx/D. The results are shown in Fig. [51 In 
zero transverse field the ground-state is degenerate and 
its energetics is governed by the Ising interaction V^^. 
The application of a small transverse field Bx lifts the 
degeneracy, with the splitting between the ground-state 
and the first excited state corresponding to the state with 
spins flipped. In that regime the most important in- 
teraction remains V^'' and the gap 5E is found to be 
proportional to Bx/D (inset of Fig. [5]), as suggested by 
the arguments leading to Eqs. ^ and By turning 

on Bx/D to larger values, the transverse field eventually 
becomes stronger than the dipolar interactions. At that 
point, the perturbative low Bx regime^ is no longer valid 
and the gap 5E is no longer proportional to Bx. How- 
ever, Fig. [5] shows that, even for high transverse fields. 
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FIG. 2: (Color Online) Comparison between the S — 1 and 
ScH ~ 1/2 models for a given sample (e.g. realization of 
disorder) of A'' = 9 spins: gap 5E / {D\/N) as a function of 
the transverse field Bx/D. 



we observe a good agreement between the S" = 1 and the 
effective S'ofr = \ description. 

Interestingly, for a specific realization of disorder, in 
Fig. [21 we note a local maximum in 5E around B^/D k, 
0.0025, followed by a local minimum, before 5E starts di- 
verging with increasing ■ We investigated the origin of 
this behavior and found that it can be understood as aris- 
ing from the B^ dependence of Czz ex (1 — ^(iJj^/rio)^) vs 
Cxo oc Bx/^o, both for small B^/D. Obviously, if this is 
the case, the random distribution of the magnetic ions in 
the sample must play a crucial role in the position (and 
even the existence) of this local maximum/minimum fea- 
ture. The structure of 5E vs B^ is controlled by the 
C^v parameters, but not only: there is also a prcfactor 
coming from the dipolar interaction which is proportional 
to rf^rf,. If one takes an extreme case in which all the 
magnetic ions are aligned on a line along the z direc- 
tion, the resultant interaction is 0, and there is no dip 
in the curve. To confirm this scenario we show in Fig. [3] 
6E/{D^/N) as a function of the transverse field B^ for 
twenty different disorder configurations for = 6. One 
sees that the majority of curves do not show these local 
maximum/minimum features and, as shown by the inset 
of Fig. [21 the average of SE over those twenty realizations 
of disorder reveal no such max/min structure. 

Having demonstrated the one-to-one correspondence 
between the 5=1 and the effective SeS = \ model 
for various (specific) realizations of disorer, we now pro- 
ceed to check the scaling with system size for (|<^i?|) pre- 
dicted by Eq. (|4]) for the 5 = 1 model and also check 
that it it agrees with the one for the effective SoS = \ 
model The results for both models are shown in Fig. [J 
The average gap (l^-El) was computed over 1000 samples 
which, for each system size of N spins, we rcnormalizc 
as {\5E\)/{D^fN), and plot for both models {S = 1 and 
SeS = \) as a function of the transverse field B^- As 



2.5e-06 



2e-06 



S 1.5e-06 

z 

Q 

S le-06 



5e-07 



_4xl0" 



3x10" 



2x10' 



1x10" 



Q- -©Average /> 



0.001 0.002 0.003 0.004 

B /D : 




0.001 0.002 0.003 0.004 0.005 
B /D 



FIG. 3: (Color Online) Random variations of the disorder 
configurations for a N=6-spin system for twenty realizations 
of disorder. Gap 5E/{Dy/N) as a function of the trans- 
verse field Bx/D. Depending on the disorder configuration 
the curves exhibit a local maximum and a local minimum. 
The thin (black) curve in the main panel shows the mini- 
mum/maximum structure of SE vs Bx for a specific realiza- 
tion of disorder. This structure disappear after taking the 
average as shown by the thick (red) curve joining the filled 
(red) circles. The monotonous behavior for the average of 5E, 
already for 20 samples, is emphasized in the inset. 



showed in Ref. [25| there exist a regime for which the 
spin 5 = 1 model obeys ^^^^ (x scahng. Indeed, 
for the 5=1 case (closed symbols), we clearly observe 
in Fig. [H a good collapse of the curves for the various 
system sizes with this linear behavior. 

One can see that at higher B^ , the scaling relation for 
different system size N , as well as the proportionality of 
the gap (|(5£'|) with Bx starts to break down. As ex- 
plained above in the context of Fig. [21 this comes from 
the fact that the transverse field term in the Hamilto- 
nian is larger than the dipolar interaction V.iY ■ Thus 
the droplet picture is not valid and neither arc the scal- 
ing nor the proportionality relations in Eq. ([5]) fulfilled. 
In Fig. [31 we also show the results for the effective 
5oft = \ model (open symbols), demonstrating the agree- 
ment with the results for the 5 = 1 model, even when the 
{SE/\/N cx Bx/D) regime breaks down. This confirms 
the correctness of the conclusion based on Eq. (|13p . and 
that SE is the same for both the 5 = 1 and the SeS = 1/2 
models. 



V. CONCLUSION 

We have shown how to rigorously derive an effective 
spin— i Hamiltonian to describe the problem of induced 
random fields in a spin glass model with strong single-ion 
Ising anisotropy and subject to a transverse magnetic 
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tive ^eff — I Hamiltonian. 



field. We discussed the relation of this problem with 
that of the LiHoa;Yi_3^F4 material in a magnetic field 
transverse to the Ho'^+ Ising spins^. We have shown, 
both analytically and numerically, that the use of such 
a model give results in full quantitative agreement with 
previously reported perturbation theory calculations on a 
"large" spin S model with strong anisotropy .^^'^^i^^ How- 
ever, the large hyperfine interactions present in the real 
LiH0j;Yi_2:F4, and which have been ignored here, must 
ultimately be considered in order to obtain a good quan- 
titative understanding of the low-temperature regimci^ 

The approach of Refs. [25l[27l | proceeds via the 
Rayleigh-Schrodinger perturbation theory, the one in 
Ref. [26| and presented in Section IIIII above relies on 
the effective Hamiltonian approach. To low order in the 
quantum 7i_L term, the two approaches have been shown 
to give identical results. However, the emergence of in- 
duced random fields is much more apparent in the spin-i 
effective model approach. The C^^ coefficients needed to 
construct the effective Hamiltonian are easily calculated, 
providing an ability to investigate the evolution of 5E and 
^ with Bx beyond the linear term and to arbitrarily high 
order in B^. Such high order perturbation theory would 
be more cumbersome to construct when proceeding via 
a direct Rayleigh-Schrodinger perturbation scheme. The 
crucial step connecting the perturbation theory method 
and the elfectivc S'off = \ approach is in the determi- 
nation of the B^ dependence of the C^^ transformation 
parameters in Eq. ([S]). It is the neglect of this B^ de- 
pendence of the spin interactions in the Scfi = \ model 
investigated in Refs. 27ll28| that seemingly led their au- 
thors to argue for the quantitative inadequacies of the 
•S'off — 5 approach. 

We note that, in a general case where the V^^^ spin- 
spin interactions are not much smaller than _D, higher- 



order perturbation theory calculations must be carried 
out to derive an effective Hamiltonian. The physical 
result would be that virtual transitions to the excited 
states lead to an admixing of those states with the 
low-energy sector. This effect was recently discussed 
in Ref. [4a|, where it was shown that such interaction- 
induced quantum mechanical effects are seemingly neg- 
ligible for LiHoa:Yi_a:F4. Tliis makes difficult to under- 
stand the advocated phenomenon of quantum mechanical 
entanglement proposed in Ref. to explain the peculiar 
behavior of the very dilute LiHoa;Yi_a;F4 {x = 0.045). 
However, as a counter-example and for a different mag- 
netic rare-earth system, we note that it was recently 
found that such interaction-induced admixing can dra- 
matically change the low-energy physics.^ 

With the contributions of Refs. ^25|i2j|2311il3|| and 
the clarification presented herein, it may be that the be- 
havior of dilute LiHo3;Yi_a;F4 in a transverse field, both 
in the random ferromagnetic and spin glass regimes, are 
now somewhat understood. This impression would seem 
to be further corroborated by recent experimental stud- 
ies which provide evidence for the manifestation of in- 
duced random fields for LiHo2;Yi_a;F4 with x = 0.44 and 
Bx > 0^ Yet, there are many questions still opened re- 
garding the physics of this material for x < 20%: Is 
there a dipolar spin glass phase over a reasonably wide 
range of dipole moment concentration, either theoreti- 
callyii^i^ or experimentall y "'''^'^'^ ? What are the physi- 
cal objects giving rise to the peculiar coherent dynamics 
at low-temperature for samples with low Ho'^^ concen- 
tration (see Refs. j6HTol|)? Even for pure LiHoF4, what 
is the microscopic explanation for the discrepancy be- 
tween experiment and Monte Carlo simulations for the 
temperature vs transverse-field phase diagram for small 
Bx near the classical paramagnetic phase boundary?—!^ 
Are the phenomena found in zero and nonzero Bx for 
LiHo3;Yi_3;F4 also observable in other Ising systems 
which possess either Kramers or non-Kramers rare-earth 
magnetic ions;^ or where the hyperfine interactions, so 
important for Ho^+ ionSf^i^iiS^ may be much less signif- 
icant3? While it is interesting that LiHo2;Yi_a;F4 in a 
transverse field becomes a rare, if not the first physical 
realization of a random-field Ising model in a ferromag- 
netic setting^ '^^1^^'^^ it would seem that this is a small 
part of the challenges offered by this material, with ap- 
parently more left to understand than has so far been 
understood. 
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